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Abstract 

An effective, low-energy, field theory of s-wave quarkonia, constituent 
heavy quarks, and gluons is constructed which is manifestly gauge invari- 
ant. The interaction Lagrangian has the form of a twist expansion, as 
typically encountered in hard processes, and involves derivatives of arbi- 
trary order. The parameters in the interaction are related with the non- 
relativistic wave function, and standard results for QQ inclusive decays 
and radiative transitions are shown to be easily recovered. The light-cone 
gluon momentum distribution at very small x is calculated and shown to 
be uniquely determined by the non-relativistic wave function. The distri- 
bution has a part which goes as x _1 logx, ie. is more singular than the 
usually assumed 1/x behaviour. The fragmentation function for a virtual 
gluon to inclusively decay into an rj c or % is also calculated. We find that 
the emission of low momentum gluons makes this process quite sensitive 
to assumptions about the binding of heavy quarks in quarkonia. 

1 Introduction 

Heavy quarkonium is traditionally modeled as a non-relativistic, colour singlet 
bound state of a Q Q pair with a static coulomb potential at short relative dis- 
tances, and some confining type of potential at long distances. The binding en- 
ergy, which is small relative to the heavy quark masses, is taken as justification for 
low quark relative velocity as well as the neglect of explicit gluon degrees of free- 
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dom in the hadron's wavefunction. In hard processes, the entire non-perturbative 
QCD physics is buried into a single parameter: the wave function at the origin 
for the case of S-wave quarkonia, or its derivative for the case of P-waves. This 
model has been widely used in calculating the decay rates of quarkonium states, 
as well as their production in e + e _ collisions, deep inelastic processes, Z° decays, 
etc [1,2]. 

The problem with this traditional approach to quarkonium modeling is 
that gauge invariance - which is obviously fundamental to QCD because it is 
its legitimising principle - is not respected as an exact symmetry. Under a local 
gauge transformation q(x, t) — > U (x, t)q(x, t), the state normally used to describe 
quarkonia, 

J d 3 x 1 d 3 x 2 f (x i - x 2 ) q (f i, t)T q (x 2 , t) | > (1) 

does not remain invariant. In the above equation T is a space-time independent 
matrix in spin, colour, and flavour indices and f(x) is the relative wave function. 
However, one can construct gauge invariant states in the following manner [3]: 
define, 

Q(x,t) =U- 1 (x,t)q(x,t). (2) 
U(x,t) = Pe ia K d *- A ®> t \ (3) 
Then, the state constructed from Q and Q, 

J d 3 Xl d 3 x 2 f {x ! - x 2 ) Q (x 1; t) T Q (x 2 , t)\0> (4) 
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is indeed invariant for arbitrary f(x). The gluon field N 1 transforms in the usual 

way, 

A* -> U^A^U - -U'^U. (5) 
9 

Unfortunately the operators Q(x, t) are not pure (current) quark fields; they also 
involve arbitrary numbers of soft gluons which are responsible for transporting 
colour between the quarks and for quarkonium binding. These constituent quarks 
are clearly extremely complicated objects. Therefore, to make a gauge invariant 
model of quarkonium requires more than that suggested by Eq.(l). Fortunately, 
traditional quarkonium models are not too widely off the mark because quarkonia 
are fairly small and the path-ordered integral in Eq.(4) is possibly negligible. For 
p- waves one expects that the problem would be more acute than for s-waves since 
a centrifugal barrier serves to keep the quarks apart, leading to a larger meson. 
However, to our knowledge, the validity of using a non-gauge invariant state for 
either s or p systems has not been investigated. 

2 The Model 

Heavy quarks can be considered as external sources placed in a gluonic vacuum 
which undergoes non-perturbative fluctuations, and results in a modification of 
the potential type of interaction between quarks in quarkonium [4]. For large 
quark mass m, one can consider a Q Q pair localized at a relative distance R such 
that m^$> 1/R and hence the relative momentum of quarks p <C m. This together 
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with the assumption that the quarkonium system is weakly bound, ensures that 
the non-relativistic approximation is valid. At the same time, we would like 
perturbative methods to be applicable, i.e. a s {p) <C 1, and hence that R <C Aq C d ' 
We shall assume that charmonium systems (bottomonium is obviously better) 
fulfil the requirement of being sufficiently heavy, yet also sufficiently small and 
weakly bound. 

We would like to construct an effective theory of quarkonia, constituent 
quarks, and gluons which respects gauge invariance. Towards this end, consider 
an elementary pseudoscalar meson field (f>(x), representing an r\ c or rjb meson for 
example, which interacts with quark fields Q(x) according to Q^^Qcf). This is 
clearly wrong since quarkonia are extended, weakly bound, systems which can be 
formed only when the heavy quark and antiquark happen to have small relative 
velocities. In contrast, a point coupling gives an amplitude for meson formation 
independent of relative velocity. To remedy this situation, and introduce the ap- 
propriate non-locality, consider an effective interaction with an arbitrary number 
of derivatives: 



Here a n are dimensionless numbers, to be determined later, and M is the quarko- 
nium mass. Before we show that this leads to conventional formulae for decays 
etc, we ask what is its gauge-invariant generalization. This is too unwieldy in 
general. But as discussed earlier, perturbation theory holds in this model and 
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so we can meaningfully consider a gauge-invariant model at the one gluon level. 
The general form of this must be: 

i D ■% D 



C-PS = ~ i Yl a nQl. 
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n=0 

oo 



M 2 



iD-iD \\J 



E K + iQ^ ^tt^ Q<j>. (6) 



, * Ulw \\~lP~) j 



sym 



where D is the usual covariant derivative, 



and the symmetrization braces are defined by 

{ }sym= (iD-idT F^+UD-lDY ( / [) ■ i /) ) : ! 7''"' ( / D ■ i D 



The interaction in Eqs. 6-8, consisting of an infinite tower of operators grouped 
together by symmetry, is, in a sense, a twist expansion of the type encountered 
in hard processes. Here the "hard momentum" is the quarkonium mass. Eq.6 
is complete at the leading twist level; other terms added on to it will be sub- 
dominant. The model leads in a straight-forward manner to Feynman vertices^. 
These are discussed below. 



■'"Eq.^) contains derivatives and therefore Hps 7^ — Lps. The presence of derivatives leads 
to additional terms upon quantization. To illustrate, suppose C ~ gip^^ipd^Q. Then the 
hamiltonian contains a term proportional to g 2 (%jj^%f) 2 . Quartic and higher self-couplings can 
be neglected at the order of accuracy of our calculations. 
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Al) PS-Quark Vertex (Fig. la). The vertex factor for pseudoscalar 
coupling to quarks is 

l5 F(p 2 ) (9) 

where, 

^)=z4y , (io) 

n=0 X" 1 / 

and ]f — \ (j>2 — ViY is the relative momentum. 

Do we have any intuition about F(p 2 ) ? Since (pi+p 2 ) 2 = M 2 , it follows that 
Ap 2 = 1p\\-1p\— M 2 approaches zero for e = 2m—M approaching zero, i.e. the weak 
binding limit. It is therefore reasonable to expect that F{p 2 ) is steeply peaked 
around p 2 = 0. These explanations are confirmed in the next section, where it will 
be shown that F{p 2 ) can be expressed directly in terms of the non-relativistic 
wavefunction of the quarks. 

A2) PS-Quark-Electric Gluon Vertex (Fig. lb). The vertex factor 
for coupling to an E\ gluon originates from expanding out the covariant deriva- 
tives in the first term of Eq.(6) and keeping a single gluon operator only. Taking 
the matrix element indicated in Fig. lb and organizing the term suitably leads to 
a rather nice and compact form, 

^75£- pG(p,q). (11) 
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The function Q (p, q) is most easily expressed in terms of the dimensionless vari- 
ables y and z, 



F(y) - F(z) 



(12) 



y - z 



where, 



y = 




m 2 



z = 




(13) 



m 2 



Note that this vertex is directly expressible in terms of F{p 2 ) and hence, as 
we shall see, in terms of non-relativistic quark wave function. This is a simple 
and direct consequence of gauge invariance, and involves no new assumptions. 
This is not true for the other gluon vertex discussed below. 

A3) PS-Quark-Magnetic Gluon Vertex (Fig.lc). This vertex 
follows from systematically expanding the second term in Eq.(6), keeping only 
one gluon, and taking the matrix element indicated in Fig.lc. This yields, 



. 9 

77T 



75 C/ii/ 



(14) 



with, 



F(y) - F(z) 



(15) 



y - z 



and, 




(16) 
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F(p 2 ) is in principle different from F(p 2 ) although one can expect a similar func- 
tional dependence. 

The symmetry of operators will be different in different mesons. Since we 
shall deal with J/ ip decays, it is useful to consider the generalization of the 
pseudoscalar results. All vertices in Fig.l are immediately applicable to the 1~~ 
system by substituting 75 — > —i^ and contracting with the meson polarization 
vector. In the limit of large M, the spin-spin interaction is weak and therefore 

F PS (p 2 )=Fy(p 2 ). 



3 Conventional Limit 



Now that the Feynman vertices for the model have been made explicit, several 
calculations can be done straightforwardly. But first, to understand the physics 
of F(p 2 ), consider the lowest order diagram (Fig. 2a) contributing to the electric 
form-factor. The contact term (Fig. 2b) involves a higher power of p 2 and is 
therefore neglected. Imagine that only the quark has electric charge and the 
antiquark is uncharged. The amplitude, to leading order in the quark relative 
momentum p and for small photon momentum g M , is 



A=2eM 2 e*.P 



* d A p 






F(p 2 )F 










2 

—m 2 +ie 


(P-\P) 2 - 


-m 2 +ie 


(p+\P-qj -m 2 +ie 



(17) 
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The p° integral may be performed by keeping only the contributions of poles in 
the vicinity of small p°. This yields, for q 2 = 0, 

F{p 2 ) 



A= -2iee*.P J d 3 p 



(18) 



(2tt) 3 /2mV2 ( e + p2/ m ) 
On the other hand, if we consider a charged spinless particle scattering from an 
e.m. field, this has an amplitude equal to —2iee*.P. This allows us to identify the 

factor in the brackets in Eq.(18) with the non-relativistic quark wave function, 

F(?) 



(2tt) 3 / 2 M 1 /2 ( e + p7m) 



nip) 



(19) 



As a consistency check, we calculate the rj — > 27 decay in the model defined 
by Eq.(6). In the cm frame both photons have large energy M/2 and therefore 
the contact term in Fig. 3b, which is sharply damped by the form factor Q, does 



not contribute. Fig. 3a and its crossed version yield for the amplitude: 



where, 



1- 



I 



d 4 p 



A = 2e 2 Me a ^ & e* a {q 1 )e*Jq 2 )q ll q 2S X. 



F(p 2 ) 



(20) 



-m 2 +ie {p — \P) -m 2 +ie {p+\q2~\qi) -m 2 +ie 



[21) 



Performing the p° integration as before, using Eq.(19) to relate F{p 2 ) to the 

wave function, doing the final phase-space integration, and summing over colours, 

one arrives at the standard expression for r\ — * 27 decays, 

1 R(0) | 2 



r^-»2 7 = 12a 2 - 



M 2 



f22) 
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Only R(0) enters the expression, which is natural enough since the two quarks 
annihilate only when very close together (and a check of the model in Eq.(6)). 
Gauge invariance of the quarkonium state does not play a significant role in this 
process. 

The emission of a soft gluon or photon, as in the Ml transition J/ip — > 
+ 7 (Fig. 4), does bring forth the issue of gauge invariance in an important way 
because the contact diagrams (Figs. 4b, 4c) are unsuppressed. The amplitude for 
the process is, 

A = -At^e a ^ s e* a (q)e^P)P,q s (l dir + Icon) ■ (23) 

Here M, P M , e^(P) refer to the J/ip. The "direct", or conventional term, follows 
from (Fig. 4a) and its crossed version, 



X dir = iM- 



(2tt)< 



(p+|-P) — m 2 +ie {p — \P) —m 2 +ie (p+^P — qj —m 2 +ie 

(24) 



Performing the p° integration and keeping only the contributions from the poles 
in the lower half plane near p° m yields, 



(2-) ;! (e + f/m) 2 



r j£p 
J (2tt) 

= J drr 2 e^Rj^{r)R n {r). (25) 



In the above, we have kept only the leading order term in the photon energy q°. 
In the dipole approximation, the exponential factor is unity and if the two mesons 
had identical wave functions then one would simply have X dir = 1. Because of 
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the hyperfine splitting this deviates from unity, and a typical (model-dependent) 
value f| is X dir = 0.987. Another set of model parameters used by the same 
authors yields 0.984, 0.920. From the amplitude, Eq.(23), the decay width is 
readily seen to be, 



|2 



3 \m, 

« 2.41 ke V. (26) 

This differs substantially from the measured width, 1.11 ± 0. 35keV. This is 
a well-known problem with the usual charmonium model and a wide range of 
explanations exist for the factor of 2-3 discrepancy. These include relativistic 
corrections, mixing effects, anomalous quark magnetic moment, etc. References 
to these may be found in a recent review by Schuler ||. 

The contact terms in Figs. 4b and 4c, which are required by gauge invariance, 
may also be calculated quite straightforwardly, 

d A p G(p,q)F(p) 



X = Ai 



(27T) 



(p + \P — |g) - m 2 + is 



(27) 



(p - \P + \qj - m 2 + ie 
Since this is a correction term, it is adequate to take Fj/^ = F v = F and 
Gj/ip — Q-q — Q ■ The p° integration gives, keeping only the contribution of nearly 
coincident poles near the origin, 

d 3 p Q{p,q)F(p) 



Icon M2J (27r)3 £+p 2/ m 



2V2 

'vrM 3 / 2 



J dpp 2 g(p,q)R(p). (2* 
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Since |<f|<^|p|, it is adequate to replace Q in Eq.(15) by, 

~..m 2 dF .. 

At this point one needs to confront the issue: what is F equal to ? We have seen 
that F is related to the electric charge distribution and is expressible in terms 
of the n.r wave function (Eq.(19)). It is possible to show that F is related to 
the magnetic response of the system and can, in principle, also be found from a 
NRQM calculation. But this is not immediately useful as these calculations have 
been done only for static quantities. Instead, we make the physically plausible 
assumption that F(p) = £F(p) where £ is a scale factor. Substituting this into 
Eqs.(29) and (28) yields, 

| roc . \ d 

Icon = - -£M ^ dp (e + p 2 /m) R(p)—pR(p) 



Any given quarkonium model allows for the calculation of e and < 4^ >• For 
definiteness, assume a gaussian wavefunction of the type exp(— p 2 /2/3 2 ), which 
yields < 1/p 2 >= 2/ ' 1 . Fitting to the i] c decay rate gives (3 2 w 0.413 Gev 2 . With 
m c = 1.65 Gev, a commonly used value for the charm quark mass, it follows that 
Icon = — 0.65£. Now suppose that the entire (large!) discrepancy between the 
measured decay width and the conventionally calculated value can be attributed 
to the magnetization term which, as we have argued, symmetry requirements 
force us to include in the Lagrangian, Eq.(6). The total decay rate is, 

T (^Si -> 7 + 1 So) = ^«eg (j-^j q (X dir + X con f 
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= 1.11 ±0.35 (forJ/^j^j + r] c ). (31) 

This suggests that £ is a number around 0.33 — 0.66. With the magnetization term 
thus determined, one can make physical predictions for quarkonium processes 
involving soft photons or gluons. An application to gluon fragmentation into 
heavy quarkonium will be described in the section after next. 

4 Gluonic Distribution 

In the present model it is possible to calculate the light-cone distribution of low 
momentum gluons in a weakly bound heavy quarkonium system. It will be shown 
that requiring gauge invariance of the hadronic state implies the existence of a 
term which goes as x _1 logx, which is more singular than the x~ l dependence 
calculated by Brodsky and Schmidt [0] using simple perturbative arguments for 
positronium. An explicit expression for the coefficient of the x~ x log a; term can 
be provided in terms of the non-relativistic wave function. Although calculable, 
we shall not worry about the x~ l terms as this is anyway a theoretical exercise - 
stable quarkonium targets unfortunately do not exist and so gluonic distribution 
are not directly measurable. 

The starting point |5| is the formula for the gluon momentum distribution 
inside a spinless hadronic target, written as a correlation of operators on the light 
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cone, 

G{x) = x / — e lXx < P | ^(0)^ (Xn) \ P> . (32) 
Here k^ and n M are two null vectors k 2 = n 2 = with k.n = 1, k~ =n + = 0, and 

Only the transverse components of the gluon field are involved. Inserting a 
complete set of states between the two operators, and limiting the outgoing X to 
two quarks only, yields, 

G(x) = x J [dp 1 }[dp 2 \8(x - q ■ n) \< p x p 2 \ A* | P >| 2 . (33) 

A summation on physical gluon polarizations (i = 1, 2), as well as colour indices 
(a = 1, 8), is implicit. The measure [dp] is, 

We now concentrate upon calculating the matrix element in Eq.(33). If we 

limit our interest to the terms most singular in x, it turns out that the term 

corresponding to radiation from a quark line (Fig. 5b) can be ignored^. The 

emission of a magnetic gluon from the vertex (Fig. 5c) is also subdominant - 

this follows because the amplitude, Eq.(14), vanishes as q —>■ 0. The dominant 

contribution comes from the emission of an electric gluon. Keeping just this term, 
2 This is the diagram which yields the x~ l behaviour in the work of Brodsky and Schmidt 

0. 
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we have from Eq.(ll), that, 



< PlP2 | A" | P >= 2^-u{p 1 ) lb v{p 2 )^p v Q. (35) 

The light-cone propagator is used in the above, 

D»» = -<T + ^ + Q — (36) 
q ■ n 

Keeping only the most singular term at small x, as well the lowest order term in 
the quark relative momentum p^ = \ (pi — p 2 ), and performing the polarization 
sum, the squared matrix element is calculated to be, 

D<»l^|f>p4¥« ! (») 

Again for small x, 

g A = _ X 2 M 2 _ g 2 ? (3g) 

[dpi] [dpi] S(x-q - n) = ^ p+^Y^ > ( 39 ) 

where K± = p 1± + p 2 ±. This yields, 

n \ - 16 9 2 1 f d P +d2 P^ K i- K lP +2 n2 (An , 
U[X) ~ xP+ J (2tt)« ( X 2 M 2 + K\r [ } 

From Eqs. (12-13) Q is related to the wave function eq (19) and its derivative 
through 

„ dF 

(41) 

where t is a dimensionless variable, 

(p 2 + \q 



,2 



2 



K\ p\ 2p +2 \ 

= -x 2 -rf. (42) 
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Performing the angular integration yields for the gluon distribution, 



G(x) 



\ogx 




(43) 



x 



where A is a positive constant determined by the non-relativistic quark wave 
function, 



This is the main result of this section. It shows that demanding gauge invariance 
of the hadronic state has a profound effect upon the distribution of low momentum 
gluons. 



Our final application of the model developed in this paper is to calculate the 
rate of fragmentation of gluons into quarkonia. Gluon fragmentation refers to 
the process of converting highly virtual gluons into hadronic physical states. The 
calculation is done in two parts. First, the fragmentation function is calculated 
at the scale of the heavy quark mass, and, second, it is evolved perturbatively 
from low to high virtualities. If one assumes that perturbative QCD is valid 
even at the scale of the charm quark mass, then a first principles calculation of 
fragmentation becomes possible. Recently Braaten and Yuan[9] have performed 
such a calculation for gluons fragmenting into ^So and 3 Si heavy mesons. 

In the first part of this section we show that the result of the calculation 




(44) 



5 Gluon Fragmentation. 
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of Braaten and Yuan |J can be exactly replicated by considering the "direct" 
diagrams (Fig.6a) implied by the model. The only difference is that our calcu- 
lation can be performed entirely in field theoretical language, which is perhaps 
an advantage. In the second part, we show that long wavelength magnetic glu- 
ons emitted by the contact diagram (Fig. 6b) augment the previous contribution. 
We remind the reader that, in the present model, "long wavelength" neverthe- 
less means a wavelength sufficiently small for perturbative QCD to be valid: as 
discussed earlier there is a hierarchy of scales, Aq C d *C q <C p <C M. 

The starting point of the calculation is the expression for the unpolarised 
gluon fragmentation function into a specific quarkonium state with momentum 
P 1 '. The reference frame is chosen to be the rest frame of the hadron, 



The notation here is identical to that in the previous section, i.e. k 2 = n 2 = 0, etc. 
The expression for D(z) follows from duplicating the analysis of Jaffe and Ji[10] 
and replacing for quark operators with gluon operators. Since D(z) involves only 
the "good" components of the gluon field, it is a twist two quantity. A sum on 
the unobserved states X is implied. To see more clearly the physical meaning of 
Eq.(45), put | PX >= C\P) | X >, where C ] {P) creates a meson of a given 
type. Using completeness of the states | X > gives, 



This makes apparent that D(z) is essentially the probability of finding a specific 



V ' 2 J 2tt 



< | A l (0) | PX >< PX | A\Xn) | > . 



(45) 




(46) 
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hadron with + component of momentum equal to zq + . The transverse momentum 
of the incoming gluon is integrated over. 

Now consider the production of a state from the process in Fig. 6a. After 
calculating traces, the amplitude for the process is, 

< PX | A" | >= l -^A» dir , (47) 

A v dir = -%g 2 me^P^l a Z dir , (48) 
d 4 p F{p-\l) 



dir 



(27T) 



(49) 



—rn 2 +ie (p—\q} —m 2 +ie (p-\-\q— lj —m 2 +ie 
The crossed diagram doubles the above value of X. Since the form-factor, which 

is essentially the wave function, restricts the relative quark momentum to small 
values, it is apparent from the denominators in Eq.(49) that the dominant con- 
tribution to the integral comes from the region around q 2 ~ 4m 2 . Performing 
the p° integral as in the previous applications, and using Eq.(19), yields for the 
direct amplitude, 

< PX | A" | > dir = ^-^^D^s^l^Fr. (50) 

where s = q 2 is the mass of the fragmenting gluon, and we have set M 2 Am 2 . 
The sum over unobserved states in Eq.(45), which amounts to an integration over 
the gluon momentum /, leads to, 

1 g A jj 2 (0) [™ ds [(1 - 2z + 2z 2 ) s 2 - 8sm 2 z + 16m 4 ] 
{Z) ~ 6^ ~JF~ J M y z ^ ( s - M 2 f { ' 
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The lower limit of integration follows from setting the minimum value of q\ 
I 2 , = in, 



M 2 z „ 
z 1 — z 



(52) 



A colour factor of 1/12 has been included in Eq.(51). Performing the integration 
yields, 



D(z) 



c? s I R(Q) I 2 

3tt M 3 



3^-2^ 2 + 2(l-z) log(l 



(53) 



This is precisely the result of Braaten and Yuan ||, ie. Eq.(8) of their paper. 

The contact diagram of Fig. 6b is calculated similarly and it has the Lorentz 
structure given in Eq.(48) with, 



1, 



d A p 



M 2 J (2tt) 4 



V + 3l) ~ m + l£ 



(p-\q)' 



— m 1 + ie 



(54) 



Only the vertex Eq.(14-16) is involved; electric gluons do not contribute here. 
Since the magnetic form factor Q restricts p to small values, from the two denom- 
inators in the above integral it is evident that the major contribution comes from 
small values of the outgoing gluon momentum I. Performing the p° integration 



gives, 



< PX I A* 1 I >„, 



d 3 p = 

\f':s I " ' J i2;;) : ' 



(55) 



Since we have chosen our reference frame as the rest-frame of the produced meson, 
it follows that, 



l\= ( s - M 2 ) /2M 



(56) 



21 



Using Eq.(29) yields, 



d 3 p ~ . . m 2 r d 3 p 1 dF 



(2tt) 3 2 7 (2tt) 3 p dp 

1 . „ (m_ . e 



- i M3 (^ + ^ Mi/2 r rfr ^ (r) j- (57) 

Thus the total amplitude is, 



< PX | A" | > dir + < PX | A" | > c 
8e 2 



M 1 / 2 s (s - M 2 ) 



^I^e^e'^ (V(0) - ±x(0)) , (58) 



where, 



X (0)=^(0) + — / drr^r). (59) 
2 Jo 

Using Eqs. (58-59) to calculate the fragmentation function yields, 

D ^ = ^A z ){m-\s^))\ (60) 



where, 



and 



■q(z) = 3z-2z 2 + 2(l -z)log(l - z) (61) 



S(0) = v^x(O)- (62) 



The first term in Eq.(60) is the direct term and is the same as Eq.(53), while 
the second is the contact term. To numerically estimate D(z), we shall make 
the same assumption as in section 3, i.e. that F{p) = £F(p) and use the same 
gaussian wave function. D(z) can then be expressed as, 
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DW = 3^^ )jR2( ° )(1 " a)2 ' (63) 
where a varies from 0.35 to 0.71 as £ goes from 0.33 to 0.66, the range estimated 
in section 3. It is clear from Eq.(63) that except for a scale factor (1 — a) 2 the 
fragmentation function at the initial scale, as well as after evolution, is identical 
to that calculated by Braaten and Yuan |§. The scale factor, however, causes a 
substantial decrease in the magnitude of D(z). 



6 Summary 

We have presented in this paper a low-energy, effective, gauge- invariant theory 
wherein the fundamental degrees of freedom are quarkonia, quarks, and gluons. 
The interaction has the form of a twist expansion familiar from hard processes 
and consists of towers of operators grouped together according to their symmetry. 
The quarkonium mass plays the role of the "hard momentum". The arbitrary 
number of derivatives in the theory serve to bring in the appropriate amount of 
non-locality or, equivalently, a form-factor in momentum space which embodies 
the extended structure of the meson. This form factor, which is the basic input 
into the model, was shown to be directly related to the non-relativistic quark 
wavefunction, a quantity calculable in any given potential model. 

An important consequence of gauge-invariance is the emergence of Feyn- 
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man vertices representing the direct gluon- quark- ant iquark-meson interaction. 
These vertices have a substantial effect upon certain heavy-meson phenomena. 
For example, the radiative M\ transition, J /if) — > i] c + 7 receives an additional 
contribution from one such vertex. This could help explain why the usual decay 
calculations invariably overestimate the decay rate by a factor of 2-3. As another 
example, we have calculated the light-cone momentum distribution of gluons 
in heavy quarkonia. Although these distributions are probably of no practical 
interest, nevertheless the present model does have some interesting theoretical 
consequences. We find that the QQG ^So vertex, which contributes to G(x), not 
only gives the x~ l behaviour but also has a x~ l logx part which is more singular 
at small x. 

As the final application of our model we considered the fragmentation of 
gluons into ^So mesons. If we ignore the contact (gauge) diagrams, then the results 
of Braaten and Yuan || are exactly recovered. But including these diagrams 
leads to a downwards rescaling of the their results by an amount which could 
be substantial. A proper calculation depends upon knowledge of the magnetic 
formfactor, called F(p 2 ) here, which in principle could be determined from a quark 
model calculation that includes states of arbitrary excitation. Finally, we remark 
that the model discussed in this paper is extendable to p states as well. This 
would be interesting because the centrifugal barrier keeps the quarks relatively 
further apart and thus makes the issue of non gauge-invariance of the meson state 
more acute. 
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Figure Captions 



1. Vertex factors for the coupling of pseudoscalar quarkonium to quarks and 
gluons. a) Quarkonium - quarks coupling; b) and c) coupling to electric 
and magnetic gluons respectively. 



2. Contributions from (a) direct and (b) contact terms to the form factor of 
'So. 



3. Contributions from (a) direct and (b) contact terms to the r\ — > 27 decay 
rate. 



4. Contributions to 3 Si — > 7+ 1 So decay rate from: a) direct diagram; b) con- 
tact diagram with 7 emanating from 3 Si vertex; c) contact diagram with 7 
emanating from 1 Sq vertex. 



5. Probing the gluon distribution in a 1 So meson: a) gluon and unobserved fi- 
nal state hadrons X; b) gluon emission from a quark line; c) gluon emission 
from x Sq vertex. 



28 



6. Gluon fragmenting into tJc^Sq). a) Direct diagram, and, b) contact diagram 



Figures may be obtained by writing to 
hafsa%png-qau%sdnpk@sdnhq.undp.org 
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